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Abstract 

We derive the rogue wave solution of the classical massive Thirring model, that describes nonlinear optical pulse propagation in 
Bragg gratings. Combining electromagnetically induced transparency with Bragg scattering four-wave mixing, may lead to extreme 
waves at extremely low powers. 
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1. Introduction 

Extreme wave phenomena appear in a variety of scientific 
and social contexts, ranging from hydrodynamics and oceanog¬ 
raphy to geophysics, plasma physics, Bose-Einstein condensa¬ 
tion (BEC), financial markets and nonlinear optics d-llll. His¬ 
torically, the first reported manifestation of extreme or rogue 
waves is the sudden appearance in the open sea of an isolated 
giant wave, with height and steepness much larger than the av¬ 
erage values of ocean waves. A universal model for describing 
the dynamics of rogue wave generation in deep water with a fiat 
bottom is the one-dimensional nonlinear Schrodinger (NLS) 
equation in the self-focusing regime. The mechanism leading 
to the appearance of NLS rogue waves requires nonlinear inter¬ 
action and modulation instability (MI) of the continuous wave 
(CW) background Indeed, the nonlinear development of 
MI may be described by families of exact solutions such as the 
Akhmediev breathers m, which are recognized as a paradigm 
for rogue wave shaping. A special member of this solution fam¬ 
ily is the famous Peregrine soliton lO, which describes a wave 
that appears from nowhere and disappears without a trace. Ex¬ 
treme waves that may be well represented by the Peregrine soli¬ 
ton have recently been experimentally observed in optical fibers 
0, in water-wave tanks ifTOl and in plasmas ifTTIl . 

Moving beyond the one-dimensional NLS model, it is im¬ 
portant to consider extreme wave phenomena in either multi¬ 
dimensional or multicomponent nonlinear propagation. Vector 
systems are characterized by the possibility of observing a cou¬ 
pling of energy among their different degrees of freedom, which 
substantially enriches the complexity of their rogue-wave fam¬ 
ilies. Recent studies have unveiled the existence of extreme 
wave solutions in the vector NLS equation or Manakov system 
C3 ca, the three-wave resonant interaction equations na, 
the coupled Hirota equations C3 and the long-wave-short- 
wave resonance ca. 
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In this Letter, we present the rogue wave solution of the clas¬ 
sical massive Thirring model (MTM) C3, a two-component 
nonlinear wave evolution model that is completely integrable 
by means of the inverse scattering transform method 1201 - 1221 . 
The classical MTM is a particular case of the coupled mode 
equations (CMEs) that describe pulse propagation in periodic 
or Bragg nonlinear optical media Ea-Ezl. Eurthermore, the 
CMEs also appear in other physical settings. In particular and 
relevant to rogue waves, they describe ocean waves in deep wa¬ 
ter for a periodic bottom 1281 . As such, the search for novel 
solution forms of these equations including rogue waves, pro¬ 
vides understanding of nonlinear phenomena and leads to appli¬ 
cations beyond optical systems. In this respect, benefiting from 
the result |25l|29l that many MTM solutions (including sin¬ 
gle and multi-solitons and cnoidal-waves) may be mapped into 
solutions of the CMEs, provides a tool used in several works 
including ocean waves |[30l, BEC ED and metamaterials (321. 

After discussing the analytical rogue wave solution in sec¬ 
tion Q, in section 0 we numerically confirm its stability, and 
show that it may also be applied to describe the generation of 
extreme events in the more general context of the CMEs. Ei- 
nally, in section Q we discuss the physical implementation of 
MTM rogue waves by using coherent effects in resonant non¬ 
linear media, such as electromagnetically induced transparency 
(LIT), which may lead to the giant enhancement of cross-phase 
modulation (XPM) with the simultaneous suppression of self¬ 
phase modulation (SPM). 

2. Analytical solution 

Let us express the MTM equations for the forward and 
backward waves with envelopes U and V, respectively, as 

Ue = -ivV--\V\^U 
y 

y = -iyu--iuiA. (1) 

y 
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Here the light-cone coordinates 77 are related to the space 
coordinate z and time variable t by the relations = dt cd^ 
and dr^ = df-cd^, where c > 0 is the linear group velocity. Even 
though the arbitrary real parameter v can be rescaled to unity, 
we find it convenient to keep it for dimensional reasons. 

The rogue waves travel over the following CW background 

Uo = ae‘‘>’ , Vo = -be‘‘^ (2) 


where, with no loss of generality, the constant amplitudes a and 
b are real, and the common phase 0 (^, rj) is 

ct> = +pr], a = bi-, 13 = . (3) 

a V by 

Up to this point we consider the two amplitudes a,b free 
background parameters. It can be proved that rogue wave so¬ 
lutions of Eqs. O exist if and only if the two amplitudes a, b 
satisfy the inequality 

0 < ab < . (4) 


By applying the Darboux method to the MTM 13^ . one obtains 
the following rogue wave solution 




) , V = -Z7e'^4(l-4*—) (5) 


AiQ 2 . 


with the following definitions 


qi = Oi(liq)q02 , q2 = 02(1 - iq)qOi , 

q = ^q + bx*^- ( 6 ) 

and 

ju = \qi\^+ \q2\^+ (i/p)i\qif, 

In the expression ([^, which is the analog of the Peregrine solu¬ 
tion of the focusing NLS equation, the free parameters are the 
two real background parameters a, b, which are however con¬ 
strained by the condition 0 , and the two complex parameters 
61 , 62 , while the parameter;^ is given by the expression 




b^ 


-(1 + ip) 


y 

a(l - ip) 


( 8 ) 


Expression 0 of the rogue wave solution may be simplified by 
fixing the reference frame of the space-time coordinates. The 
general solution ^ may then be obtained by applying to this 
particular solution a Lorentz transformation. 

According to the last remark above, we now provide the 
rogue wave solution in terms of the space z = c(^ - q) and time 
t = (^ -\-q) coordinates directly. By rewriting the CW phase ^ 
in these coordinates, one obtains 


0 

OJ 


7 7^/1 ab^b 

kz-ojt, k= —(1 - —)(- - -) , 
2 c y^ a b 

y ab a b 

-;z(l- ^)(t + -) , 

2 y^ b a 


(9) 


where k is the wave number of the background CW. Setting 
a = b means choosing the special frame of reference such that 
^ = 0. Note that the other possibility a = -b does not satisfy 
the condition that p is real (see (|7])). Erom a physical standpoint, 
the CW background solution with a = b corresponds to a non¬ 
linear wave whose frequency oj = -y(l - a^/y^) enters deeper 
inside the (linear) forbidden band-gap 0 / < as its intensity 
grows larger. A linear stability analysis of the CW background 
solution Q shows that it is modulationally unstable for pertur¬ 
bations with a wavenumber (for details, see El). 

Note that modulation instability gain extends all the way to ar¬ 
bitrarily long-scale perturbations (albeit with a vanishing gain), 
a condition which has been refereed to as ’’baseband instabil¬ 
ity”, and that is closely linked with the existence condition of 
rogue waves in different nonlinear wave systems (e.g., the Man- 
akov system, see Ref. lfTSl ). It is also interesting to point out 
that, outside the range of existence of the rogue wave solution 
that is for > y^, the background is unstable with respect 
to CW perturbation with a finite (nonzero) gain (see Ref. El)- 
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Figure 1: Evolution in time and space of the intensities in the forward and 
backward components of the rogue wave solution 


By using translation invariance to eliminate the parameters 
61,62 , one finally ends up with the following expression of the 
MTM rogue wave solution 


U = ae — [1- q*(q-^i)] , 

p 

V = — [I - -q*(q-\- /)] 

p* p 


( 10 ) 


where 


OJ = 

P = 
P = 


-y(l - -v) , q = 


— [ip(z - zo) - c(t - to)] 
yc 


— - 1 


2 M 2 + (l+ 2 Im^)(l--), 
P 


( 11 ) 


where zo and to are arbitrary space and time shifts, respec¬ 
tively. Note that a further simplification may come from rescal¬ 
ing z, t, U, y by using the length scale factor S = -ycja^. 

In Eigj^we show the dependence on space and time of the in¬ 
tensities |f/p and |yp of the forward and backward rogue waves 
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( p^ . Here we have set v = -1, c = I, a = 0.9, tQ = 2 and 
Zo = 3.5. As can be seen, the initial spatial modulation at r = 0 
evolves into an isolated peak with a maximum intensity of about 
nine times larger than the CW background intensity. The corre¬ 
sponding contour plot of these intensities is shown in Fig|^ 



Figure 2: Contour plot of the intensities of the forward and backward waves as 
in Fig[T] 


3. Numerical results 

In order to verify the spatio-temporal stability of the rogue 
wave solution ( p^ over a finite spatial domain, we numerically 
solved Eqs. Q with y = -1, c = 1, and using the initial (i.e., at 
f = 0) and boundary (i.e., at z = 0 and z = L) conditions given 
by the exact expression ( p^ , with the same solution parameters 
as in Figs[T]|^ and with L = 5 ^ 4S. 

Fig|^ displays the numerically computed intensities of the 
forward and backward waves: as it can be seen, there is an 
excellent agreement with the analytical solution. This confirms 
the stability and observability of the rogue wave solution 
in spite of in the presence of the competing background MI. 
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Figure 3: Numerical solution corresponding to the analytical solution in Figs[^ 

At this point it is quite natural, and interesting, to numerically 
check whether the initial conditions of the rogue wave solution 
may induce also the generation of a rogue wave as modeled 


by the CMEs for pulse propagation in nonlinear Bragg gratings 
1^ or ocean waves with periodic bottom |[28]| . Indeed, the 
CMEs 

= -ivV-^(\Vf + cr\Uf)u, 

= -ivU-^{\U\^ + cr\V\^)v. (12) 


differ from the MTM Eqs. Q by an additional SPM term, its 
relative strength being expressed by the coefficient cr. We have 
numerically solved Eqs.lp^ with initial and boundary condi¬ 
tions as given by expression (1^). As shown by Fig|^and Fig|^ 
quite surprisingly even in the case of a comparatively large SPM 
contribution, one still observes the generation of an extreme 
peak with nearly the same intensity and time width as in the 
MTM case. However now the peak no longer disappears, but 
it leaves behind traces in the form of dispersive waves, owing 
to the non-integrability nature of the CMEs. In Figwe set 
cr = 0.3, while qualitatively very similar results are also ob¬ 
tained for larger values of cr (e.g., for cr = 0.5, which corre¬ 
sponds to the case of nonlinear fiber gratings). 
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Figure 4: Numerical solution of CMEs with cr = 0.3 and initial and boundary 
conditions as given by the analytical solution of Figs[^ 

On the other hand, in Figwe show the case with cr = -0.5, 
that is relevant to water wave propagation in oceans with a pe¬ 
riodic bottom ||28]| . Interestingly, in the hydrodynamic case the 
formation of a first rogue peak where both components have 
nearly the same amplitude as in the MTM case Fig|^ is fol¬ 
lowed by wave breaking into secondary, and yet intense peaks 
in each of the two generated pulses that travel in opposite direc¬ 
tions. 


4. Discussion and conclusions 

For the strict applicability of the analytical MTM rogue wave 
solution ( p^ of section Q, it is necessary that SPM can be ne¬ 
glected with respect to XPM. This situation occurs whenever 
the forward and backward wave envelopes have different car¬ 
rier frequencies, say, oju and ojy, and their frequency difference 
Aoj = oJu - ojy is close to a resonant frequency of the nonlin¬ 
ear medium Ea. Whenever the forward and backward waves 


3 



















5. Acknowledgments 



Space z Space z 


Figure 5: As in Figs[^ with cr = -0.5 . 


are coherently coupled in a four-level atomic system by a CW 
pump field, a giant enhancement of the strength of XPM is also 
possible via the BIT effect l^ - ITTl . with no competing SPM. 
In addition, BIT brings the important benefit of zero linear ab¬ 
sorption losses, and it is only limited by the residual two-photon 
absorption, which can be neglected for relatively short interac¬ 
tion distances. 

Bor maximizing the interaction length of the counter- 
propagating waves, diffraction effects can be suppressed by us¬ 
ing an atomic vapor cell containing an hollow-core photonic 
band gap fiber (PBGB). A scheme of Bragg soliton generation 
in a coherent medium exhibiting BIT has been previously dis¬ 
cussed, however counter-propagating pulses at the same probe 
frequency were considered (381 . A dynamic refractive index 
grating for coupling forward and backward signal waves at dif¬ 
ferent frequencies may be induced by the standing wave gener¬ 
ated from the beating of two counter-propagating pump waves, 
with a frequency difference equal to A6 g>. This process is known 
as Bragg scattering four-wave mixing (BS-BWM) (391, and it 
enables noise-free frequency translation. Indeed, BS-BWM in 
the co-propagation geometry and at microwatt pump power lev¬ 
els has been recently demonstrated in Rubidium vapor, confined 
to a few centimeters long PBGB (40l . 

Note that the solution ( p^ may also apply to describe ex¬ 
treme wave emergence in the co-propagation of two signals 
at different carrier frequencies (411 . by simply interchanging 
the role of time and space variables in Bqs. Q- In this case, a 
dynamic grating may be induced by the beating of two pump 
waves in the orthogonal polarization (42l . 

In short summary, we obtained the rogue wave solution of 
the classical MTM, which extends the Peregrine soliton solu¬ 
tion of the NTS equation to the case of wave propagation in 
a periodic nonlinear medium. An implementation is proposed 
using coherent resonant wave mixing: we envisage that Bragg 
solitons and rogue waves may ultimately be observed at sub¬ 
milliwatt pumping levels by using chip-scale waveguides that 
are evanescently coupled to atomic vapors (H. 
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